Charged vortices in superfluid systems with pairing of spatially separated carriers 
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It is shown that in a magnetic field the vortices in superfluid electron-hole systems carry a real 
electrical charge. The charge value depends on the relation between the magnetic length Ib and the 
Bohr radiuses of electrons a% and holes a B . In double layer systems at filling factors v e — Vh = v 
and for a B ,a B ^> £b the vortex charge is equal to the universal value ve . 
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It is generally believed, that the vortices in supercon- 
ductors are connected with an applied magnetic field, 
while the magnetic field does not have any influence on 
the properties of the vortices in electrically neutral su- 
perfluid systems. The aim of this letter is to show that 
in superfluid systems subjected by a magnetic field the 
vortices will have a real electrical charge (the compensat- 
ing charge of the opposite sign will appear on the surface 
of the system) . In general case the charge of the vortices 
takes a fractional value. For the first time the fractional 
charge of the vortices was predicted by Laughlin [|l| for 
the two-dimensional electron gas in a quantized magnetic 
field. Then, it was established in Ref. ^ |^|, that in 
double layer electron systems with the half-filling of the 
lowest Landau levels in each layer the vortices should 
carry the charge equals to ±e/2 (here and below e is the 
absolute value of the electron charge). 

We will show that in any superfluid system the mag- 
netic field results in an appearance of the vortex charge 
proportional to the polarizability of the particles and in- 
verse proportional to their effective mass. The estimates 
show that for reachable values of magnetic fields the vor- 
tex charge is unobservable small one in superfluid phases 
of He isotopes and in Bose gases of alkali metals, while 
it can be of order of the electron charge in superfluid 
systems with pairing of spatially separated electrons and 
holes. 

The authors of Ref. 0] call the possibility of the 
electron-hole pair superfluidity in question based on the 
fact that the interband transitions fix the phase of the or- 
der parameter and result in a transition into a dielectric 
state. But it was established in Refs. || |(| that the inter- 
diction on the electron-hole superfluidity can be removed 
in systems, where the spatially separated electrons and 
holes are coupled. In these systems the interband transi- 
tions coincide with the intcrlaycr ones and usually they 
are exponentially small. The superfluid state of the pairs 
with spatially separated components has both the super- 
fluid and superconducting features. The superfuid flow 
of such electron-hole pairs is accompanied with real su- 
percurrents flowing in opposite directions. Therefore, we 
will call these systems the condenser superconductors. 

In Refs. §, | the pairing of a conducting band elec- 



tron from the one layer with a valence band hole from the 
other layer was considered. Then in a number of theoret- 
ical papers [|[ |[ 0, & ^, [l(| it was shown a possibility of 
superfluidity of pairs composed from spatially separated 
electrons and holes belonging to the conducting band. 
This possibility is realized in double layer electron sys- 
tems in a magnetic field normal to the layers for the case 
of the total filling factor vt = ^1 + ^2 = 1< During almost 
10 years there were many efforts to observe the condenser 
superconductivity experimentally |yl [l2| [H| Q . Now it 
seems that these effort have been crowned with success 

©Hi- 

The principal result can be obtained from general con- 
sideration which does not imply a concrete form of the 
Hamiltonian. Let us consider a superfluid system sub- 
jected with crossed electric E and magnetic B fields. Let 
£ is the density of the energy of the system, and II is 
the density of its generalized momentum. Note, that the 
values of £ and II do not include the energy and the mo- 
mentum of the external fields that polarize the system. 
We are interesting in a relation between the momentum 
density II and the dipole momentum density P. To es- 
tablish this relation it is convenient to consider the frame 
of reference, in which the electric field is equal to zero. 
The velocity of that frame relative to the lab frame is 
equal to 



E x B 



(1) 



The expression ([l]) is valid to within the linear order of 
the factor u/c. The condition u/c -C I is satisfied at 
E 2 <C B 2 . Below we assume this condition is fulfilled. 

Let us introduce the notation £ and n , the energy 
density and the momentum density in the new frame, 
correspondingly. The relation between the energies and 
the momenta in the lab frame and in the new one can be 
obtained from the transformations 



n = n, 



pu, 



£ = £o(n ) 



IIqu 



pu 



(2) 



(3) 
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Here p is the mass density. Expressing the momentum 
IIo in Eq. O) in terms of the momentum II we obtain 



£ = £ (n- P u) + nu-^- 



Then, we take into account that 



d£_ 



5E 



where v is the velocity of the superfiuid system. 
Eqs. (§) and 

we find the required relation 



n 



-P x B. 



(4) 



(5) 
Using 



(G) 



In general case the kinematic momentum pv is of order 
of the generalized momentum II. But, as it is shown be- 
low, for the condenser superconductors in a strong mag- 
netic field the kinematic momentum is much smaller than 
the generalized one. Consequently, in a strong magnetic 
field the term pv in Eq. (^) can be omitted. Taking into 
account, that in the superfiuid system II = nh\7<p, where 
n is the density of the particles (we consider the temper- 
ature T = 0), and if is the phase of the order parameter, 
we obtain from Eq. (|^) the following relation 



nfiVip 



-P x B. 

c 



(7) 



Let the system, lying in the (x, y) plane, is subjected 
with a uniform magnetic field, directed along the z axis. 
Taking the curl of the both parts of Eq. ((?]) , we obtain 



nh cur^V^ = — div2P, 
c 



(8) 



where div 2 P is the two-dimensional divergence, of the 
vector P. This quantity taken with its sign changed is 
equal to the polarization charge density p po \. On the 
other hand, the left hand side of Eq. (||) is nonzero only 
in case, when the quantized vortices exist in the system. 
Then 



curl 2 V</5 = 2ir S(r — r^)r 



(9) 



where nt = ±1 the upper (low) sign corresponds to the 
vortices rotating in the counter-clockwise (clockwise) di- 
rection, and the summation is over the vortex centers. 
Integrating the both parts of Eq. (|J) over an arbitrary 
area, for vortices of the same sign one finds 



± 2-xhnN,, 



H 

c 



PpoidS 



H 



(10) 



Here N v is the number of the vortices in the area S, and 
Q v is their charge. It follows from this, that vortex charge 
is equal to 



Qv , „ he 
— = ±2ir — n 
N v B 



±2iret B n. 



(11) 



In the case of the electron-hole pairing in the lowest Lan- 
dau level, the density of the pairs n is related with the 
filling factors v e = = v by the formula v — 2iri 2 B n. 
Thus, in strong magnetic fields in the superfiuid phase 
the quantized vortices carry an electrical charge q = ±ve. 

To clarify how general is the result obtained let us an- 
alyze the behavior of the condenser superconductor in a 
strong magnetic field perpendicular to the layers at small 
filling factors v, when the electron-hole pair gas can be 
considered as a weakly interacting Bose gas. Let us con- 
sider two 2D conducting layers separated by the dielectric 
layer of the width d with the electron carriers in one layer 
and the hole carriers in the adjacent layer. We consider 
rrih S> m e . In a strong magnetic field the large differ- 
ence of the masses m e and rrih may result in a situation, 
when the electron Bohr radius a e B = eh 2 /m e e 2 is much 
larger then the magnetic length £ B = (ch/eB) 1 ^ 2 , and 
the hole Bohr radius a B = eh 2 /irihe 2 can be smaller or 
larger then £ B . In this case the energy spectrum of the 
bounded electron-hole pair, which is formed due to the 
Coulomb interaction between spatially separated carri- 
ers, was found in |Q. The part of the pair energy de- 
pending on the momentum of the pair tt and the velocity 
u is equal to 



AS 



2M* 



Mb 
M* 1 



IM B 2 
mhU . 

2 M* 



(12) 



Here A/* = Mb + rrih, the pair effective mass, M B , the 
"magnetic mass" 



M R = 



4 eh 2 



'2tt e 2 £ B 

Introducing the pair polarizability 



2tt £b 



a(B) = M B ^ = 4^=4, 



one can rewrite the correction AS in the form 

2 



1 / . ,ExB 



(13) 



(14) 



(15) 



Analogous expression was obtained in 1 18 for an electri- 



cally neutral atom in crossed fields for the case of small 
magnetic fields. In that case Eq. ( |l5| ) contains the zero 
magnetic field polarizability of the atom a(0) instead of 
a(B) and the mass of the atom M instead of the mass 
of the pair M* . 

Replacing the momentum tt with the operator —iKV, 
from Eq. ( p"5| ) we obtain the Hamiltonian of the electron- 
hole pairs. In the low density limit, when the size of 
the pair is much smaller then the distance between the 
pairs and the exchange effects are inessential, the pairs 
can be considered as true bosons. At low temperatures 
the rarefied Bose gas should form a superfiuid state. The 



3 



superfluid phase can be described by the order parameter 
*. The order parameter satisfies the equation 



ih- 



9* 
dt 



1 



2AL 



-iHV + a(B) 



E x B 



^p-E 2 ^ + 1 \^\ 2 -$>. (16) 



The last term in the r.h.s. of Eq.([L6|) describes the in- 
teraction between the pairs. One can show that in the 
the limit d <C Ib the interaction constant is equal to 
7 = (7r/2) 3 / 2 e 2 d 2 /e£s. The vanishing of the interaction 
constant at d = is the consequence of the exact compen- 
sation of the Coulomb forces between the pairs (compare 
with JT^I). Presenting the order parameter in the form 
vj/ = |vj/| e ^( r ), we obtain from Eq.(]l|) the velocity of the 
superfluid component 



1 



hVtp + a{B)- 



E x B 

he 



(17) 



To obtain the dipole momentum of the unit area P 
we take into account that the r.h.s. of Eq. (16) is the 
variational derivative over <3>*(r) of the energy functional 
of the Gunzburg-Landau type for superconductors. The 
derivative of that functional over the electric field E taken 
with the opposite sing is P: 



a(B) 



1 



M*c 



V</2 x B 



1*1 



The expression ( |18| ) can be rewritten in the form 
P = a(£)(E+ -v s x BWl 2 . 



(18) 



(19) 



This result means that not only the electric field E, 
but also the Lorentz force polarizes the pair, acting in 
opposite directions on the positive and negative charges 
of the pair. However, such a simple expression for P is 
valid within the linear accuracy in v s /c. 

One can note that the expressions ( |l7| ) and ([l8]) are 
in agreement with Eq. (Q) . Actually, if we exclude from 
Eqs. ( |l7| ) and (|l^) the electric field E and take into 
account that in the case considered p = to/, J*! 2 , we arrive 
to the relation (||). 

The polarization charge p po i is found by taking the 
two-dimensional divergence of the both sides of Eq. ( 18 ) . 
To calculate the derivatives in the r.h.s. of Eq. ( |l8| ) we 
take into account that the relation (|^) for the curl z V<y5 
takes place. Beside this we omit the terms containing 
the quantity V|*|. We are not interested here in the 
structure of the vortex core (we consider it as the math- 
ematical point), and out of the core the terms containing 
V|^| are small at small fields E. It allows to replace 1*1 2 
with the pair density n. Finally, we obtain 



Ppoz(r) = -a(B)n 



Mb\ 
M* ) 



div 2 E- 



(20) 



Thus, the polarization charge in a superfluid system in 
a magnetic field can be caused by the polarization of the 
medium by the electric field with a nonzero divergence, or 
by the existence of the quantized vortices in the system. 
It also follows from Eq. ((2^), that the charge of the z-th 
vortex is equal to the coefficient of S(r — r'), namely 



hB , . 
q = ±2n—a(B)n. 
M*c 



(21) 



Eq. ( pl| ) yields the vortex charge for the electron-hole 
double layer systems in an arbitrary magnetic fields. The 
same result is valid for the electron-electron double layer 
system with the substitution M* = Mb + 2m e . 

In weak magnetic fields (is ^> a|) the polarizability 
a = j(a e B ) 3 , where 7 ~ 1, and the effective mass Af* ~ 
to/, + to 6 ~ to/,. Then, the vortex charge is equal to 



± 



2-7T7 



ft B a B ne. 



(22) 



In high magnetic fields 
for a(B) one obtains 



> e B > a h B ) using Eq. 



M 



M F 



q = ±2nl B n e = ± — ■ — ve 



M, 



M, 



(23) 



Finally, in ultra high fields (a B » £b) the effective mass 
M* = M B {l+y/2nt B /±a B ) -> M B and the vortex charge 
obeys the universal relation q = ±ve. 

Let us present here some estimates. In the magnetic 
field B = 10 T the magnetic length £ B « 80 A. In GaAs 
heterostuctures with the dielectric constant e — 13 the 
magnetic mass is very small (Mb ~ 10~ 28 g). Due to 
the smallness of this quantity the vortex charge can be 
of order the electron charge only in superfluid systems 
with the boson mass of order the electron mass Too- For 
the electron-hole double layers at m h = 0.4 m and m e = 
0.067too we obtain M^/M* w 0.2 and the vortex charge 
is equal to q ss 0.2i/e. For the electron-electron double 
layer system with the same m e we find Ms/M* w 0.5 
and the vortex charge q w 0.5ve The derivation of q 
from the universal value is connected with that the ratio 
£3/0,% ~ 0.8 and the strong inequality is <C a B does not 
satisfy in this case. 

The number of the vortices and their spatial distribu- 
tion depends on the electric field E. Integrating the both 
sides of Eq. p^) along a certain contour, we find 



M* <pv s dl = 2nhN v 



^B^Exdl. (24) 



Here N v is the number of the vortices inside the contour 
(we consider the vortices of the same vorticity). It fol- 
lows from Eq.(p4[), that the velocity v s can be reduced 



4 



under the appearance of the vortices - in such a way we 
lower the kinetic energy of the system. When the vortex 
distribution is considered as a continuous one, the vortex 
density n v (r) can be introduced 



N v -- 

Putting the r.h.s. of Eq.(| 
(pj) and (H) the relation 



n v (r)dr. 



(25) 



to zero we find from Eqs. 



n v (r) 



a(B)B 
2i:hc 



div 2 E. 



(26) 



In the case of a condenser superconductor of a disk 
shape and for the electric field directed along the disk 
radius and independent of the angle (p, it is reasonable 
to assume, that the vortices situated on the circumfer- 
ences centered at the disk center. If the radius of the 
j-th circumference is pi and the number of the vortices 
on it is Ni, then it follows from the definition of n v (r) 



n v (r)rdrd9. 



(27) 



This equation determines the relation between the 
quantities pi and Ni. The values of pi can be found 
from the assumption that the average distance between 
the vortices on a given circumference (which is equal 
2irpi/Ni) coincides with the average distance pi+i - pi 
between the vortices belonging to the adjacent circum- 
ferences. It yields 



Ni = 



2npi 



Pi+i 



(28) 



Eqs. ([26|)-(|28|) allow to find the values of Ni and pi up 
to a factor of order of unity. 

A macroscopic number of the vortices with equal vor- 
ticites can also emerge in the absence of the electric field. 
It is realized when besides the uniform field B z there is 
an extra field B T with div2B T ^ (B T is parallel to the 
plane of the structure). Indeed, one can show (compare 
with ) that in the field B T the energy of the pair of 
spatially separated electron and hole is equal to 



£ = 



1 , ed „ n9 
(tt -\ z x B r ). 

2M„ V c ' 



One can see that energy (g9|) differs from the expressions 
( |l5| ) only by that the induced dipole momentum a(B)~E 
is replaced with the spontaneous momentum edz. There- 
fore, the dipole momentum of the unit area can be ob- 
tained from Eq. ( |l9|) replacing the induced momentum 
with the spontaneous one. Then, taking the divergence 
of P , we find 



Ppoi{r) 



a{B)Bn 



— div 2 B r 

c 



2 2irhni8(v — r^) 

(:■!(» 



It follows from this expression that in this case the vortex 
charge is equal the value found above and in the continu- 
ous limit the vortex density is n v (r) = (ed/27r?ic)div2B r . 



At nonzero temperatures the charged vortices will 
emerge in condenser superconductors in a fluctuation 
way, in similarity with the same phenomena in a thin He- 
ll film. The circumstance that the vortices are charged 
does not influence in the first approximation on the ther- 
modynamic features of the system. It is connected with 
that the Coulomb correction to interaction between the 
vortices falls down much faster (by the power law) then 
bare logarithmic interaction between them. The last one, 
as is well known, results in a Kosterlitz-Thouless transi- 
tion. Since the sign of the vortex charge is in one to one 
correspondence with the sign of the vorticity, at tem- 
peratures below the Kosterlits-Thouless temperature the 
vortex-antivortex pairs should be electrically neutral. At 
temperature above the Kosterlits-Thouless the vortices 
and antivortices decouple, and free electrical charges ap- 
pear. It reveals itself in a principal change of the conduct- 
ing properties of the system under the phase transition. 
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